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The fast Monte Carlo procedure of hadron generation developed in our previous work is extended 
to describe noncentral collisions of nuclei. We consider different possibilities to introduce appropri- 
ate asymmetry of the freeze-out hyper-surface and flow velocity profile. For comparison with other 
models and experimental data we demonstrate the results based on the standard parametrizations 
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of the hadron freeze-out hyper-surface and flow velocity profile assuming either a common chemi- 
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I. INTRODUCTION 

In the preceding work [l| we have developed a Monte Carlo (MC) simulation procedure, 
and the corresponding C++ code allowing for a fast but realistic description of multiple 
hadron production in central relativistic heavy ion collisions. A high generation speed and 
easy control through input parameters make our MC generator code particularly useful for 
detector studies. The generator code is quite flexible and allows the user to add other 
scenarios and freeze-out surface parametrizations as well as additional hadron species in 
a simple manner. We have compared the BNL Relativistic Heavy Ion Collider (RHIC) 
experimental data on central Au+Au collisions with our MC generation results obtained 
within the single freeze-out scenario with Bjorken-like and Hubble-like freeze-out surface 
parametrizations. Although simplified, such a scenario nevertheless allowed for a reasonable 
description of particle spectra and femtoscopic momentum correlations. This description 
can be farther improved by introducing finite emission duration and extending the table 
of the included resonances; the single freeze-out scenario is however less successful in the 
description of the data on elliptic flow (see section III). 

The particle densities at the chemical freeze-out stage are too high (see, e.g., [2j) to 
consider particles as free streaming and to associate this stage with the thermal freeze-out 
one. In this work we have implemented as an option more sophisticated scenario of thermal 
freeze-out: the system expands hydrodynamically with frozen chemical composition, cools 
down and finally decays at some thermal freeze-out hypersurface. The RHIC experimental 
data are compared with our MC generation results obtained within this thermal freeze-out 
scenario. We do not consider here a more complex freeze-out scenario taking into account 
continuous particle emission (see, e.g., Q]). 

In present paper, we also extend the fast Monte Carlo procedure of hadron generation 
developed in our previous work [l| to describe noncentral collisions of nuclei. One of the 
most spectacular features of the RHIC data is large elliptic flow Q]. The development of 
a strong flow is well described by the hydrodynamic models and requires short time scale 
and large pressure gradients, attributed to strongly interacting systems. However, results 
of hydrodynamic models significantly disagree with the data on femtoscopic momentum 
correlations (compare 5] with, e.g., [6|), related with the space-time characteristics of the 
system at freeze-out. Usually, the hadronic cascade models underestimate the momentum 
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.g. q y, 3). 



anisotropy and overestimate the source sizes (e.g. |7|, |8|, |9|]). Some sophisticated hybrid 
models (e.g. AMPT 10j) reproduce the elliptic flow and the correlation radii but with 
different sets of model parameters. 

Successful attempts to describe simultaneously the momentum-space measurements and 
the freeze-out coordinate-space data were done in several models which make experimen- 



tal data fitting within some parametrizations of 



reeze-out hypersurface: "Kiev-Nantes" 



model jjj, "Blast- Wave" parametrizations [ll, 12, 13], "Buda-Lund" hydro approach JjJ 
All these approaches use the hydro-inspired parametrizations of freeze-out hypersurface and 
help in understanding the full freeze-out scenario at RHIC. 

In this article we analyze the RHIC data at y/s NN = 200 GeV and try to use the same set 
of the model parameters for the description of both the momentum-space observables, i.e. 
transverse mass (m t ) spectra and p r dependence of elliptic flow, and freeze-out coordinate- 
space observables, i.e. /^-dependence and azimuthal angle ($) dependence of the correlation 
radii. The chemical composition of the fireball was fixed in our previous article [l| by the 
particle ratios analysis. 

The paper is organized as follows. Section [III is devoted to the description of main 
modifications of the model [l] needed to take into consideration noncentral collisions. In 
section II I II the example calculations are compared with the RHIC experimental data. We 
summarize and conclude in section IIVI 

II. FREEZE-OUT SURFACE PARAMETRIZATIONS 

The extension of our MC generator to noncentral collisions demands mainly the mod- 
ifications of freeze-out hypersurface parametrizations (Sec. V of Ref. [lj]) and does not 
practically influence the generation procedure itself (Sec. VI of Ref. [1]). Therefore we focus 
on these modifications only considering the popular Bjorken-like and Hubble-like freeze-out 



parametrizations respectively used in so-called blast wave [11( and Cracow [15( models as 
the example options in our MC generator. Similar parametrizations have been used in the 
hadron generator THERMINATOR Q. 

As usual, in the Bjorken-like parametrization, we substitute the Cartesian coordinates t, 
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z by the Bjorken ones 



r = {t 2 -z 2 f' 2 , ^ = iln^, (1) 



and introduce the the radial vector r = {x, y} = {r cos0, rsin0}, i.e., 

= {r coshr], r, r sinh r/} = {r cosh rj, r cos 0, r sin 0, rsinh^}. (2) 

For a freeze-out hypersurface represented by the equation r = r(r], r, 0), the hypersurface 
element in terms of the coordinates n, r, becomes 

,o dx a dx l3 dx 1 , , , , N 

rf ^ = e ^^dT dT] d ^ (3) 

where e MQ!( g 7 is the completely antisymmetric Levy-Civita tensor in four dimensions with 
e 0123 = — eoi23 = 1- Generally, the freeze-out hypersurface is represented by a set of equations 
r = rj (rj, r, 0) and Eq. (jHJ) should be substituted by the sum of the corresponding hypersur- 
face elements. For the simplest and frequently used freeze-out hypersurface r = const, one 
has 

g? 3 <7 m = n^o = t d 2 rdi] {cosh rj, 0, 0, — sinhry}, 

d 3 a = Td 2 rdr], (4) 
?2 M = {coshr/, 0, 0, sinhr/}. 
In noncentral collisions the shape of the emission region in the transverse (x-y) plane can 
be approximated by an ellipse (as usual, the z-x plane coincides with the reaction plane) 



The e 
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lipse radii R x (b) and Ry(b) at a given impact parameter b are usually parametrized 11 



19| in terms of the spatial anisotropy e(b) = (R 2 — i? 2 ,)/ (R 2 X + R 2 ) and the scale factor 
R s (b) = [(R 2 x + R 2 )/2Y/ 2 , 



R x (b) = R 8 (b)y/l-e{b), R y {b) = R 8 (b)y/l + e{b). (5) 
Then from the ellipse equation, 

follows the explicit dependence of the fireball transverse radius R(b, 0) on the azimuthal 
angle 0: 



'I - e 2 (b) 

R{bA)=R s {b)^ y ; (7) 
a/1 + e[b) cos 20 
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particularly, R(b, 0) = R x {b) and R(b,n/2) = R v {b). To reduce the number of free parame- 



ters, we assume here a simple scaling option 



20) 



R 8 (6) = R s {b = 0) v/l-e s (6), (8) 

where R s (b = 0) = i? is the fireball freeze-out transverse radius in central collisions. It means 
that the dimensionless ratio R s (b) / R s (0) at the freeze-out moment depends on the collision 
energy, the radius Ra of the colliding (identical) nuclei and the impact parameter b through 
a dimensionless e s (b) only. It should be noted that both e s (b) and the fireball freeze-out 
eccentricity e(b) are determined by the eccentricity e (b) = 6/(2i2^) of the elliptical overlap 
of the colliding nuclei at the initial moment, when 

Rs(b) __ # s (&)imtial 



Rs(b = 0) 



'l-co(6). (9) 



e(6)=6o(6) R A 

Since e s (0) = e(0) = eo(0) = 0, one can can assume that e s (b) ~ e(b) at sufficiently small 
values of the impact parameter b. It appears that the use of the simple ansatz e s (b) = e(b) 
allows one to achieve the absolute normalization of particle spectra correct within ~ 10% 
up to b ~ Ra (see section ITlI CI) . 

If the system evolution were driven by the pressure gradients, the expansion would be 
stronger in the direction of the short ellipse x-axis (in the reaction plane), where the pressure 
gradient is larger than in the direction of the long ellipse y-axis (see, e.g., 6|). The typical 
hydrodynamic evolution scenario is shown in Fig. [TJ During the evolution, the initial system 
coordinate anisotropy eo(b) is transformed into the momentum anisotropy 6(b). According 
to the hydrodynamical calculations, the spatial eccentricity almost disappears and the mo- 
mentum anisotropy saturates at rather early evolution stage before freeze-out. As we do not 
trace the evolution here, we will consider the spatial and momentum anisotropies e(b) and 
5(b) as free parameters. 

For central collisions the fluid flow four-velocity u^(t,x) = j(t, x){l, v(t, x)} = 
7(i, ^){1, v±(t, x), v z (t, x)} at a point x and time t was parametrized |l] in terms of the 
longitudinal (z) and transverse (_L) fluid flow rapidities 

1 l + v z (t,x) 1 1 + v±(t,x) cosh r] u (t,x) 

2 l — v z (t,x) 2 1 — v±(t, x) coshr] u (t, x) 

where v± = \v±\ is the magnitude of the transverse component of the flow three- velocity 
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v = {v±cos(j) u ,v±sin(f) u ,v z }, i.e., 

u^(t,x) = {cosh p u cosh Tj u , sinh p u cos U , sinh p u sin U , cosh p u sinh 
= {(1 + u^) 1 / 2 cosh^, (1 + w 2 ) 1 / 2 sinhr^}, 

w± = 7'5j_ = r y±coshr] u v±, r y± = coshp u . However, unlike the transverse isotropic 
parametrization (0 U = 0), now the azimuthal angle n of the fluid velocity vector is not 



necessarily identica 



parameter 5(6) 
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;o the spatial azimuthal angle 0, because of the nonzero flow anisotropy 



19] 



u^(t,x) = {7^ cosh p u cosh r] u , a/1 + 5(6) sinh p u cos 1 
a/1 - 5(b) sinh p n sin 0, 7^ cosh p u sinh 7/ u } , 



where 



(12) 



70 = y 1 + 5(6) tanh p n cos20, (13) 



, . - 5(6) , > 

tan0„ = \ j - - — - tan0. (14) 

The transverse flow rapidity p u is related to p u by: 

u± = sinhp u = a/1 + 5(b) cos 20 sinhp u . (15) 

Note, that for 5(6) = (i.e. n = 0), Eq. (TT2"j) reduces to Eq. (ITTj) which was applied in 



Refs. 
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2l| . In Ref. 19|, 5(6) is obtained by fitting the model prediction to the measured 
elliptic flow coefficient v%. 

Further we assume the longitudinal boost invariance Vu — V: which is a good ap- 
proximation for the highest RHIC energies at the midrapidity region. To account for the 
violation of the boost invariance, we have also included in the code an option corresponding 
to the substitution of the uniform distribution of the space-time longitudinal rapidity 77 in 
the interval [— r) max , ?7 max ] by a Gaussian distribution exp(—r] 2 /2Ar] 2 ) with a width param- 
eter At] = r] max . The presence of the "oscillation term" yl + 5(6) cos 20 in the transverse 
component u± of the flow velocity in Eq. ffl5|) allows us to use the simple linear profile for 
p u without introduction of the additional parameters for each centrality (6) unlike other 
models, namely: 

^=im pT * {b = 0) ' (16) 

where p™ ax (6 = 0) is the maximal transverse flow rapidity for central collisions. At such 
normalization and 5(6) > e(6) the maximal transverse flow (u±,p u ) is achieved at = 0, 



i.e. alone; it should be according to the hydrodynamic scenario described above 

(Fig. [1]). (although p u has a maximum at — 7r /2!) 

Here one should note that the "popular parametrization" of transverse flow rapidity used 
in Ref. [11] (and implemented as an option in our MC generator also): 



Pu = r[p (b) +p 2 {b) cos: 



(17) 



where 



'r cos ( 



r 



r sini 



Rai J V Ry ) R(b, > 



(18) 



is the "normalized elliptical radius", po(b) and p 2 {b) are the two fitting parameters, is close 
to our parametrization and gives the similar results for observables under consideration. In 
parametrization of Ref. Uj the boost is perpendicular to the elliptical subshell on which 
the source element is found: tan0 M = (R^/R 2 ) tan0 = (1 — e)/(l + e)tan0 and 5(b) = 
2e(6)/(l + e 2 (&))- It is interesting to note that for sufficiently weak transverse flows, p u < 1, 
considered here, one can put sinh p u ~ p u and obtain our parametrization from that of 
Ref. [11] by substitutions 



p (b) p™ x (b = 0) 



I f ^cos20 u 



l + 5(b) cos 20. 



(19) 



R(b,cf>) R s (b) Po (b) 

Thus, in the case of moderate transverse flows, one can obtain the same result either by 
fixing the direction of the flow velocity vector but allowing for the azimuthal dependence 
of the flow rapidity or by allowing for arbitrary direction of the flow velocity vector but 
assuming azimuthally independent flow rapidity. 

At t = const, the total effective volume for particle production in the case of noncentral 
collisions becomes 

»27T 



where (n^u 11 ) 



d ap(t, £)« M (£, x) = t 

a(t,x) JO 



cosh p u \/l + 6(b) tanh 2 p u cos 20 . 



R(b,<p) fflmn 

(n jU u M )rdr / dn, 



(20) 



We also consider the Cracow model scenario 15] corresponding to the Hubble-like freeze 



out hypersurface th = (t 2 ~ x 2 — y 2 — z 2 ) 1 ^ 2 = const. Introducing the longitudinal space-time 
rapidity n according to Eq. (Op) and the transverse space-time rapidity p = sinh -1 (t/th), 



one has 



22] 



x 



T#{cosh n cosh p, sinh p cos 0, sinh p sin 0, sinh n cosh p}, 



(21) 



th = t~b I cosh p. Representing the freeze-out hypersurface by the equation th = 
Th(t), p,<f>) = const, one finds from Eq. (j3J): 

d 3 a = rfj sinh p cosh pdr/dpdcj) = r H dr]d 2 r, 

(22) 

n^it, x) = x^(t,x)/T H . 



With the additional flow anisotropy parameter 5(b) the flow four- velocity is parametrized 

as 

//'■ 



111 11 

(t, x) = {7? cosh p cosh T], y/l + 5(b) sinh p cos 1 



a/1 — 5(b) sinh p sin , 7^ cosh p sinh 77 } , 

where 



(23) 



7 ^ = ^1 + 5(b) tanh 2 p cos 20. (24) 
The effective volume corresponding to r = r H sinh p < R(b, 0) and 7/ min < 77 < 77 max is 

/• r2lT rR(b,4>) fljmax 

V eff = / d 3 a lx (t,x)u' 1 (t,x) = t h / d0 / (n^)rdr \ dt] (25) 



fr(t,x) JO JO Jrj: 



m 1 n 



with 

(n^u 11 ) = cosh 2 p( + tanh 2 p cos 20 

, v v x (26) 

- tanh 2 p(y/l + 5(b) cos 2 + y/l - 5(b) sin 2 0)J ~ 1 + o(5 2 (6)). 

Our MC procedure to generate the freeze-out hadron multiplicities, four-momenta and 
four-coordinates for central collisions has been described in detail in Ref. 1)]. For noncentral 
collisions, only the generation of the transverse radius r is slightly different, taking place in 
the azimuthally dependent interval [0, R(b, <f>)]. 



III. INPUT PARAMETERS AND EXAMPLE CALCULATIONS 



A. Model input parameters 

First, we summarize the input parameters which control the execution of our MC hadron 
generator in the case of Bjorken-like and Hubble-like parametrizations, and should be spec- 
ified for different energies, ion beams and event centralities. 

1. Thermodynamic parameters at chemical freeze-out: temperature T ch and chemical 
potentials per a unit charge pb, Ps, Pq- As an option, an additional parameter 7 S < 1 
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takes into account the strangeness suppression according to the partially equilibrated 



distribution 
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Up* ; T, in, ls ) = 5 ^ , (27) 

T^expCK-^/^il' 

where nf is the number of strange quarks and antiquarks in a hadron i, p*° is the 
hadron energy in the fluid element rest frame, g»j = 2Jj + 1 is the spin degeneracy factor 
Optionally, the parameter 7 S can be fixed using its phenomenological dependence on 



the temperature and baryon chemical potential 25|. 



2. Volume parameters: the fireball transverse radius R(b = 0) (determined in central 
collisions; in noncentral collisions we use the scaling option (!H|9|) to recalculate R(b) 
from R{b = 0)), the freeze-out proper time r and its standard deviation At (emission 



duration) 26] 



3. Maximal transverse flow rapidity p™ ax (& = 0) for Bjorken-like parametrization in cen- 
tral collisions. 

4. Maximal space-time longitudinal rapidity r/ max which determines the rapidity interval 
[ — iJmiB! )]max] in the collision center-of-mass system. To account for the violation of 
the boost invariance, we have included in the code an option corresponding to the 
substitution of the uniform distribution of the space-time longitudinal rapidity rj in 
the interval [— ?7 max , f/ max ] by a Gaussian distribution exp(—r] 2 /2Ar] 2 ) with a width 



parameter Ai] = r\ max (see, e.g., |20L |27|) 



5. Impact parameter range: minimal b m \ n and maximal 6 max impact parameters. 

6. Flow anisotropy parameter 5(b) in Bjorken-like and Hubble-like parametrizations (or 
Po(b) and p2(b) in the "Blast- Wave" parametrization of Ref. [ill]). 

7. Coordinate anisotropy parameter e(b). 



8. Thermal freeze-out temperature T (if single freeze-out is considered, T = T ). 

9. Effective chemical potential of ir + at thermal freeze-out fi e ^ th (0, if single freeze-out 
is considered). 

10. Parameter which enables/disables weak decays. 



TABLE I: Model parameters for central Au + Au collisions at "v/^jvtv = 200 GeV. Chemical freeze- 
out parameters: T ch =0.165 GeV, ^=0.028 GeV, /I 5 =0.007 GeV, J1 Q = - 0.001 GeV. 

T th , GeV 0.165 0.130 0.100 

7.0 7.2 8.0 

2.0 2.0 2.0 

9.0 9.5 10.0 

0.65 0.9 1.1 

0. 0.10 0.11 



TABLE II: Model parameters for Au + Au collisions at v^attv = 200 GeV at different central- 
izes. Chemical freeze-out parameters: T ch =0.165 GeV, /Ib=0.028 GeV, £ts=0.007 GeV, J1q= - 
0.001 GeV. Thermal freeze-out parameters: T th =0.1 GeV, th =0.11 GeV. Volume parameters 
determined in the central collisions: R(b = 0) =10.0 fm, r =8.0 fm/c, p™ ax (& = 0) = 1.1 



centrality 


c=0-5 % 


c=5-10 % 


c= 10-20 % 


c=20-30 % 


c =30-40 % 


c =40-60 % 


frmin/ Ra 


0. 


0.447 


0.632 


0.894 


1.095 


1.265 


bmax/RA 


0.447 


0.632 


0.894 


1.095 


1.265 


1.549 


e(b) 











0.1 


0.15 


0.15 


5(b) 


0.05 


0.08 


0.12 


0.25 


0.34 


0.36 



The parameters used to simulate central collisions are given in Table [B The param- 
eters determined in central collisions for T th =0.1 GeV: r=8.0 fm/c, R(b = 0) = 10. fm, 
Ar=2.0 fm/c; p™ ax (6 = 0) = 1.1 (3-th column in Tabled]) were used to simulate Au+Au 
collisions at \fs NN = 200 GeV at different centralities. The additional parameters needed 
only for noncentral collisions are given in Table [III 

B. Different chemical and thermal freeze-outs 

Since the assumption of a common chemical and thermal freeze-out can hardly be justified 
(see, e.g., j^), we consider here a more complicated scenario with different chemical and 
thermal freeze-outs. 



r, fm/c 
At, fm/c 
R(b = 0), fm 

pr x (b = o) 

,.eff th 
r"rr 
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The mean particle numbers iV* h at thermal freeze-out can be determined using the fol- 
lowing procedure In our preceding article [l| the temperature and chemical potentials 
at chemical freeze-out have been fixed by fitting the ratios of the numbers of (quasi) stable 
particles. The common factor, V^, and, thus, the absolute particle and resonance num- 
bers was fixed by pion multiplicities. Within the concept of chemically frozen evolution 
these numbers are assumed to be conserved except for corrections due to decay of some 
part of short-lived resonances that can be estimated from the assumed chemical to thermal 
freeze-out evolution time. Then one can calculate the mean numbers of different particles 
and resonances reaching a (common) thermal freeze-out hypersurface. At a given thermal 
freeze-out temperature T th these mean numbers can be expressed through the thermal ef- 
fective volume V£g an d the chemical potentials for each particle species /if 1 . The latter can 
no longer be expressed in the form fa = which is valid only for chemically equilibrated 
systems. For a given parametrization of the thermal freeze-out hypersurface, the thermal 
effective volume V^g (and thus all /i* h ) can be fixed with the help of pion interferometry 
data. 

In practical calculations the particle number density p° q (T, /ij) is represented in the form 
of a fast converging series [l|: 

P?{T,Pi) = ^Tj2^eM^)K 2 (^l (28) 

71 k=l 

where K 2 is the modified Bessel function of the second order, m ; and = 2J, + 1 are the 
mass and the spin degeneracy factor of particle i respectively. 

Using Eq. (I28p and the assumption of the conservation of the particle number ratios 
from the chemical to thermal freeze-out evolution time, we obtain the following ratios for 
z-particle specie to 7r + : 

The absolute values of particles densities p^ q (T th , (if 1 ) are determined by the choice of the 
free parameter of the model: effective pion chemical potential /^ ff th at the temperature of 
thermal freeze-out T th . Assuming for the other particles (heavier then pions) the Boltzmann 
approximation in Eq. ( 128]) one deduces from Eqs. ( 1281) - ( 1291) the chemical potentials of 
particles and resonances at thermal freeze-out: 

cq/ T ch ch\ ^eq/^th eff th\ 
,.th_T-ithi ( Pi V 1 ; Pi I Pit K 1 ; Pit ) \ MM 
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The correct way to determine the best set of model parameters would be achieved by 



fitting all the observables together as it was suggested in Ref. [27( , but for our MC-type model 
it is technically impossible. For the example calculations with our model at RHIC energies 
we choose T ch = 0.165 GeV and the thermal temperatures as in the analytical models which 
performed the successful fitting of RHIC data: T th = T ch = 0.165 GeV (Cracow model Jjj]) 
and T th = 0.100 GeV (Blast- Wave model III]), and some arbitrary intermediate temperature 
T th = 0.130 GeV. It is well known (see, e.g., 2j) that the pion transverse spectra at thermal 
freeze-out can be described in two regimes: low temperature and large transverse flow on 
the one hand, and higher temperature and non-relativistic transverse flow on the other hand 
(see section UlI Cp . The low temperature regime seems to be preferable because the strong 
transverse flow is expected to describe the large inverse slopes of transverse spectra of the 
. leavy hadrons (especially protons) and small correlation radii obtained at RHIC better 
3], Q • We present the calculated correlation radii in section IIIIEI 

In the considered here last version of FASTMC the new table of resonances was included. 
It contains 360 resonances and stable particles, instead of 85 ones included in the previous 
versions. This particle table is produced from the SHARE particle table excluding not 
well established resonances states. The decays of resonances are controlled by its lifetime 
1/r, there T is the width of resonance specified in the particle table, and they occur with the 
probability density Texp(— Ft) in the resonance rest frame. Then the decay products are 
boosted to the reference frame in which the freeze-out hypersurface was defined. Because 
we need to compare our calculations with data from different experiments we made possible 
to switch on/off different decays based on their lifetime (i.e. turn on/off weak decays). Only 
the two- and three-body decays are considered in our model. The branching ratios are also 



taken from the particle decay table produced from the SHARE decay table 
decays are also possible. 



28j . The cascade 



C. mj-spectra 

n 

In Fig.[2]the m r spectra measured by the STAR Collaboration [29j at — 5% centrality are 
shown for tc + , K + and p in comparison with the model calculations under the assumption 
of the common chemical and thermal freeze-out at T th = T ch = 0.165 GeV (Fig. [2£a)) and 
under the assumption that the thermal freeze-out at T th = 0.100,0.130 GeV occurs after 
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the chemical one (Fig. [2]^b, c)). 

The correction on weak decays was introduced by the STAR Collaboration in pion spectra 
. It was approximately 12% and was estimated from the measured K® and A decays. 
29| the STAR Collaboration doesn't introduce the weak decay correction in proton 



only 
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In Ref. 

spectra. To reproduce the STAR weak decay correction procedure, we excluded pions from 
and A decays from pions m^-spectra in Fig. [21 The contribution of weak decays in the 
simulated proton spectra can be estimated from Fig. [2] by comparison of the solid lines 
(protons from K® and A decays are included) and the dashed lines (without contribution 
of protons from the weak decays). The model parameters at different temperatures are 
presented in Table HI The parameters were optimized this way to obtain the good description 
of the pion m r spectra and the correlation radii. The best description of the m r spectra was 
achieved at T th = 0.100 GeV (Fig. EJ^c)). 

The same set of parameters T, p™ ax , R and r which was determined for central collisions 
(Table [TJ) was used for noncentral ones. The additional parameters of the model for non- 
central collisions were coordinate and momentum asymmetries: e and S (Table HT1). At the 
freeze-out moment we consider them as free parameters because we do not trace the evo- 
lution here. The influence of the choice of e and 5 on m r spectra averaged over azimuthal 
angle tp is negligible. The decrease of the effective volume in noncentral collisions (Eq. 120]) 
due to nonzero values of e and 5 allows us to obtain the correct absolute normalization 
of mi-spectra without introduction of the additional parameters. In Fig. [3] the mj-spectra 
measured by the STAR Collaboration 29j are shown for tt + , K + and p at centralities: 
- 5%, 5 - 10%, 10 - 20%, 20 - 30%, 30 - 40%, 40 - 50% in comparison with the model 
calculations which assume that the thermal freeze-out at T th = 0.1 GeV occurs after the 
chemical one (solid lines). It appears that the procedure described in section HT1 allows one 
to achieve the absolute normalization of pion spectra correct within ~ 13%. 



D. Elliptic flow 



30 



Following a standard procedure 
distribution in the azimuthal angle (p at mid-rapidity: 

dN dN 



311 ] we make a Fourier expansion of the hadron 



d 2 p t dy 2np t dp t dy 



1 + 2t> 2 cos2<£> + 2V4 cosA(p + ...). 



(31) 
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The elliptic flow coefficient, t>2, is defined as the second order Fourier coefficient, 

f n 27r dip cos 2(ip — ib R ) , f N , 

~ f 2n TiD d " N ' 

JO r dycUpptdpt 

where tpR is the reaction plane angle (in our generation i/jr — 0), y and p t are the rapidity 
and transverse momentum of particle under consideration, respectively. 

The value of v 2 is an important signature of the physics occurring in heavy ion collisions. 
According to the typical hydrodynamic scenario shown in Fig. [TJ the elliptic flow is generated 
mainly during the high density phase of the fireball evolution. The system driven by the 
internal pressure gradients expands more strongly in its short direction (into the direction 
of the impact parameter x in Fig. [TJ which is chosen as a "positive" direction) than in 
the perpendicular one ( "negative" direction, y in Fig. [TJ where the pressure gradients are 
smaller. Figure [TJ illustrates qualitatively that the initial spacial anisotropy of the system 
disappears during the evolution, while the momentum anisotropy grows. The developing 
of strong flow observed at RHIC requires a short time scale and large pressure gradients, 
which are characteristics of a strongly interacting system. The reason for the generation of 
i>2 at the early times is that the system should be hot and dense, when the system cools and 
become less dense the developing of the large pressure gradients becomes impossible. The 
elliptic flow coefficient, t/2, depends on the transverse momentum p t , the impact parameter 
b or centrality, as well as, the type of the considered particle. All these dependencies have 
been measured at RHIC [32]. 

The ^-dependence of v 2 measured by the STAR Collaboration [32| for charged particles 



at centralities: - 5%, 5 - 10%, 10 - 20%, 20 - 30%, 30 - 40%, 40 - 60% is shown in Fig. 1 
in comparison with our MC calculations obtained with the optimal model parameters from 
Table [III The calculations were performed under the assumption that thermal freeze-out at 
T th = 0.1 GeV occurs after the chemical one at T th = 0.165 GeV. 

The calculations under the assumption of the common chemical and thermal freeze-out 
at T th = T ch = 0.165 GeV demonstrate not so good agreement with the experimental data 
at small pt < 0.4 GeV/c for the centralities larger than 20%; irrespective of the choice of e 
and 5 one cannot get a satisfactory description in the whole p t -range (see e.g. Fig. [5]). 
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E. Correlation radii 



The parameters of the model presented in Table [J were optimized to obtain the best 
description of the pion mt-spectra and the correlation radii in the following cases: under the 
assumption of the common chemical and thermal freeze-out at T th = T ch = 0.165 GeV and 
under the assumption that the thermal freeze-out at T th = 0.100, 0.130 GeV occurs after the 
chemica. one. In Fig. the fitted correlation radii and R lmg are compared with 

those measured by the STAR Collaboration [5J. The three-dimensional correlation function 
was fitted with the standard Gaussian formula: 

CF( Pl ,p 2 ) = 1 + A exp(-i^ ut ^ ut - R 2 sidc q 2 side - flLgftong), (33) 

where q— pi —p~2 = (<?out, <?side, liong) is the relative three-momentum of two identical particles 
with four-momenta pj and pi. The form of Eq. (1331) assumes azimuthal symmetry of the 
production process [33]. Generally, e.g., in the case of the correlation analysis with respect 



to the reaction plane, all three cross terms can be significant 27|. We will consider 



this case below. We choose the longitudinal co-moving system (LCMS) as the reference 



frame 34|. In LCMS each pair is emitted transverse to the reaction axis so that the pair 
rapidity vanishes. The parameter A measures the correlation strength. 

The regime with the large temperature T th = T ch = 0.165 GeV was tested in Ref. pj. We 
have repeated this test here with the new resonances table and the additional parameter At 
(Fig. E(a), dashed line). We have found that these modifications lead to a better description 
of the correlation radii. In Fig. [6](a, bottom) (dashed line) the intercept A is larger than 
the experimental one, but taking into account the secondary pions from the weak decays 
essentially improves the description of the A (Fig. [6](a, bottom), solid line). 

In Fig. [6](b, c) we consider the lower thermal freeze-out temperatures: 0.130,0.100 GeV. 
The secondary pions coming from the weak decays were taken into account. 

It is worth to note a good description of the correlation radii (within ~ 10% accuracy) 
altogether with the absolute value of the m t spectra in the scenario with a low temperature 
thermal freeze-out of chemically frozen hadron-resonance gas. There are three important 
reasons for this success. First, a relatively small (compared with dynamic models) effective 
volume of the system ~ tR 2 that reduces the correlation radii. Second, relatively large 
transverse flow in the model that further reduces the radii. Third, rather large effective 
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pion chemical potential which is needed to describe the absolute value of the pion spectra at 
relatively small effective volumes; it reduces correlation radii at small pt and so makes their 
m t behavior flatter. This reduction happens due to vanishing of the homogeneity length of 
Bose-Einstein distribution for low-p t pions when the pion chemical potential approaches the 



pion mass (see also Ref. |35j for the analysis of the reduction of the pion correlation radii 
near the point of the Bose-Einstein condensation in static systems). We do not consider 
here the question whether such conditions could be realized in realistic dynamical models. 

It should be noted that the description of the /^-dependence of the correlation radii 
has been achieved within ~ 10% accuracy for all three considered thermal temperatures: 
T th = 0.165, 0.130, 0.100 GeV. However, at lower temperatures there is more flexibility in the 
simultaneous description of particle spectra and correlations because the effective volume 
isn't strictly fixed as it is in the case of the single freeeze-out (T th = T ch = 0.165 GeV). In 
present work, we have not attempted to fit the model parameters (T th , R, r, //® th ) since 
it is rather complicated task requiring a lot of computer time. We have performed only 
example calculations with several sets of the parameters. 

In noncentral collisions the measurement of azimuthally sensitive correlation radii pro- 
vides the additional information about the source shape. For the corresponding femtoscopy 



formalism with respect to the reaction plane see, e.g., }18l. l27j]. In the absence of azimuthal 
symmetry, the three additional cross terms contribute to the Gaussian parametrization of 
the correlation function in Eq. (|33l) : 



CF(p 1} p 2 ) — 1+A exp ( - R Q q out -R s q sidc - R\ g long - 2 R os q out q sidc - 2 _R ol g out gi ong - 2 -R sl g S id e <2Wg ) • 

(34) 

In the boost-invariant case, the transverse-longitudinal cross terms R 2 { and R^ vanish in the 
LCMS frame, while the important out-side R 2 S cross term is present. 

In the Gaussian approximation, the radii in the Eq. (1341 are related to space-time vari- 



ances via the set of equations 
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27|: 



R 2 = l/2((x 2 ) + (y 2 )) - l/2((£ 2 > - (y 2 )) cos(2$) - (xy) sin(2$), 
R 2 Q = l/2({x 2 ) + (y 2 )) + l/2((£ 2 > - (y 2 )) cos(2$) + (xy) sin(2$)) 

-2(3 ± ((tx) cos($) + (ty) sin($)) + /3 2 (P>, (35) 
R 2 = (^)-2p i m+Pf(t 2 ), 
i4 = (xy) cos(2$) - l/2((x 2 } - (y 2 )) sin(2$) + (3±((tx) sin($) - (ty) cos($)), 
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where (3i = k z /k°, (3± = k±/k° and $ = Z(kj_,b) is the azimuthal angle of the pair three- 
momentum k with respect to the reaction plane z-x determined by the longitudinal direction 
and the direction of the impact parameter vector b = (x, 0,0); the space-time coordinates 
are defined relative to the effective source center (x M ): x M = x^ — (x M ). The averages are 
taken with the source emission function S(t,x,k), 181 ]: 

i fu _ fd 4 xf{t,x)S(t,x,k) 

{f ^ X)) ~ Jd*xS(t,x,k) ' [ J 

The illustrative calculations of the correlation radii as a function of the azimuthal angle 
$ were done with the following fast MC parameters: T th = 0.1 GeV, p™ ax (b = 0) = 1.0; 
R(b = 0) =11.5 fm, r =7.5 fm/c, At =0. fm/c, e = 0.1 and 5 = 0.25. The azimuthal 
dependence of the correlation radii in different k t intervals is shown in Fig. [7J 

The R 2 S oscillates downward, in the same phase as "RHIC" source extended out of plane 



361 ]. which means the larger sideward radius viewed from the x-direction (in the reaction 



plane), than from ^/-direction (out-of plane). The source has small coordinate asymmetry 
e = 0.1, it is almost round (as in Fig. [1] step 3), however the emission zone, or "homogeneity 
region" , varies with $ because of the non-isotropic flow. 



IV. CONCLUSIONS 



We have developed a MC simulation procedure and the corresponding C++ code, that 
allows a fast realistic description of multiple hadron production both in central and noncen- 
tral relativistic heavy ion collisions. A high generation speed and an easy control through 
input parameters make our MC generator code particularly useful for detector studies. As 
options, we have implemented two freeze-out scenarios with coinciding and with different 
chemical and thermal freeze-outs. We have compared the RHIC experimental data with our 
MC generation results obtained within the single and separated freeze-out scenarios with 
Bjorken-like freeze-out surface parameterization. 

Fixing the temperatures of the chemical and thermal freeze-out at 0.165 GeV and 0.100 
GeV respectively, and, using the same set of the model parameters as for the central col- 
lisions, we have described single particle spectra at different centralities with the absolute 
normalization correct within ~ 13%. 

The comparison of the RHIC v-i measurements with our MC generation results shows 
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that the scenario with two separated freeze-outs is more favorable for the description of the 
Pt-dependence of the elliptic flow. 

The description of the /^-dependence of the correlation radii has been achieved within 
~ 10% accuracy. The experimentally observed values of the correlation strength parameter 
A has been reproduced due to the account of the weak decays. 

The analysis of the azimutal dependence of the correlation radii indicates that the source 
considered in the model oscillates downward, in the same phase as " RHIC" source extended 
out of plane. 

The achieved understanding of the reasons leading to a good simultaneous description 
of particle spectra, elliptic flow and femtoscopic correlations within the considered simple 
model could be useful for building of the complete dynamic picture of the matter evolution 
in A+A collisions. 
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FIG. 1: The typical hydro-dynamic evolution scenario. 
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FIG. 2: mt-spectra (in c 4 /GeV 2 ) measured by the STAR Collaboration [29] for tt + (circles), K + 
(squares) and p (up-triangles)at — 5% centrality in comparison with the model calculations at 
T th = 0.165(a), 0.130(6), 0.100(c) GeV , with the parameters from TablelH for protons weak decays 
are taken into account (solid lines); for protons weak decays are not taken into account (dashed 
lines). The direct tt + , K + and p contributions are shown on (c) by dotted lines. 
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FIG. 3: mf-spectra (in c 4 /GeV 2 ) measured by the STAR Collaboration [23] for tt + (circles), K + 
(squares) and p (up-triangles) at different centralities in comparison with our fast MC calculations 
at T th = 0.100 GeV (solid lines) with the parameters from Table [J and Table UU 
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FIG. 4: The ^-dependence of V2 measured by the STAR Collaboration [32|] (points) for charged 
particles at different centralities in comparison with our fast MC calculations at T th = 0.100 GeV 
(solid line) with the parameters from Table U and Table HH 
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p t (GeV/c) 

FIG. 5: The pt-dependence of V2 measured by the STAR Collaboration (points) for charged 
particles at centrality 20 — 30% in comparison with our fast MC calculations under assumption of 
the single freeze-out at T th = T ch = 0.165 GeV. The different sets of coordinate and momentum 
assymetries parameters were tried: e = 0.1, 5 = 0.3 (solid line), e = 0.2, 5 = 0.4 (dotted line), 
e = 0.1, 5 = 0.4 (solid line), e = 0.2, 5 = 0.5 (dotted line), e = 0.15, 6 = 0.4 (dashed line) 
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FIG. 6: The tt + correlation radii at mid-rapidity in central Au+Au collisions at \^s NN 



200 GeV 



from the STAR experiment jg] (open circles) and MO calculations within the Bjorkcn-likc model 
with the parameters presented in Table [J in different intervals of the pair transverse momentum kt- 
The full calculation with resonances (a) , (b) . (a) sing le freeze-out T ch = T th = 0.165 GeV, no weak 
decays (dashed line), with weak decays (solid line); (b) thermal freeze-out at T th = 0.130 GeV 
occurs after the chemical one, weak decays are taken into account (solid line); (c) the full calculation 
with resonances, weak decays are taken into account at T th = 0.100 GeV (solid line), the direct 
pions only (dotted lines). 
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FIG. 7: Simulated with FASTMC squared correlation radii versus the azimuthal angle $ of the 
7r + 7r + pair with respect to the reaction plane, 20-30 % centrality events in fcy (GeV/c) intervals: 
0.15 < kr < 0.25 (solid line), 0.25 < k T < 0.35 (dashed line), 0.35 < k T < 0.45 (dotted line), 
0.45 < kx < 0.60 (dotted-dashed line), simulation was done with the special set of parameters: 
T th = 0.1 GeV, p™ ax (b = 0) = 1.0; R(b = 0) =11.5 fm, r =7.5 fm/c, Ar =0. fm/c, e = 0.1 and 
5 = 0.25, weak decays were not taken into account. 
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